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Abstract. 

The powerful narrow jets observed to emanate from many compact accreting ob- 
jects may arise from the twisting of a magnetic field threading a differentially rotating 
accretion disk which acts to magnetically extract angular momentum and energy from 
the disk. Two main regimes have been discussed, hydromagnetic outflows, which have 
a significant mass flux and have energy and angular momentum carried by both the 
matter and the electromagnetic field and, Poynting outflows, where the mass flux is 
negligible and energy and angular momentum are carried predominantly by the elec- 
tromagnetic field. Here we consider a Keplerian disk initially threaded by a dipole-like 
magnetic field and we present solutions of the force-free Grad-Shafranov equation for 
the coronal plasma . We find solutions with Poynting jets where there is a continuous 
outflow of energy and toroidal magnetic flux from the disk into the external space. 
This behavior contradicts the commonly accepted "theorem" of Solar plasma physics 
that the motion of the footpoints of a magnetic loop structure leads to a stationary 
magnetic field configuration with zero power and flux outflows. 

In addition we discuss recent magnetohydrodynamic (MHD) simulations which es- 
tablish that quasi-stationary coUimated Poynting jets similar to our Grad-Shafranov 
solutions arise from the inner part of a disk threaded by a dipole-like magnetic field. 
At the same time wc find that there is a steady uncollimated hydromagnetic outflow 
from the outer part of the disk. The Poynting jets represent a likely model for the jets 
from active galactic nuclei, microquasars, and gamma ray burst sources. 



INTRODUCTION 

Highly-coUimated, oppositely directed jets are observed in active galaxies and 
quasars [1], and in old compact stars in binaries [2,3]. Further, well collimated 
emission line jets are seen in young stellar objects [4,5]. Different ideas and models 
have been put forward to explain astrophysical jets [6] with recent work favoring 
models where twisting of an ordered magnetic field threading an accretion disk 



acts to magnetically accelerate the jets [7]. There are two regimes: (1) the hy- 
dromagnetic regime, where energy and angular momentum are carried by both the 
electromagnetic field and the kinetic flux of matter, which is relevant to the jets 
from young stellar objects; and (2) the Poynting flux regime, where energy and 
angular from the disk are carried predominantly by the electromagnetic field, which 
is relevant to extra galactic and microquasar jets. 

Here, we discuss recent theoretical and simulation work on Poynting jets. Station- 
ary Poynting flux dominated jets have been found in axisymmetric MHD simula- 
tions of the opening of magnetic loops threading a Keplerian disk [8, 9]. Theoretical 
studies have developed models for Poynting jets from accretion disks [10-12]. 



THEORY OF POYNTING OUTFLOWS 

Consider the coronal magnetic field of a differentially rotating Keplerian accretion 
disk for a given poloidal field threading the disk. The disk is perfectly conducting 
with a very small accretion speed. Further, consider "coronal" or force-free mag- 
netic fields in the non-relativistic limit. Cylindrical (r, 0, z) coordinates arc used 
and axisymmetry is assumed Thus the magnetic field has the form B = 'Bp + B^4> , 
with Bp = BrY + B^T. . We can write B^ = -{l/r){d^! / dz), B^ = {l/r){d'^ /dr), 
where ^(r, 2;) = rA^{r,z). In the force-free limit, the magnetic energy density 
B^/ (Stt) is much larger than the kinetic or thermal energy densities; that is, the 
flow speeds are sub-Alfvenic, <^ f ^ = B^/47rp, where va is the Alfven velocity. 
In this limit, ~ J x B; therefore, J = AB [13] and consequently, 

A** = -H(t)'^ . 

^ ^ d"^ 

with A* = d'^/dr'^ — {l/r){d/dr) + d'^/dz'^, which is the Grad-Shafranov equation 
for (see e.g. [10]). 

We consider an initial value problem where the disk at t = is threaded by a 
dipole-like poloidal magnetic field. The form of H{^) is then determined by the 
differential rotation of the disk: The azimuthal twist of a given field line going 
from an inner footpoint at ri to an outer footpoint at r2 is fixed by the differential 
rotation of the disk. The field line slippage speed through the disk due to the disk's 
finite magnetic diffusivity is estimated to be negligible compared with the Keplerian 
velocity. For a given field line we have rd0/S^ = dsp/ Bp, where dsp = V dr'^ + dz^ 

is the poloidal arc length along the field fine, and Bp = B^ + B^. The total twist 
of a field line loop is 

with the sign included to give A0 > 0. For a Keplerian disk around an object of 
mass M the angular rotation is ujk — ^JCM/r^ so that the field line twist after a 
time t is 



(1) 



A0(^) = UJo t 
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where tq is the radius of the O— point, Uq = ^GM/tq, and JF is a dimensionless 
function (the quantity in the square brackets). 

The Grad-Shafranov equation (1) can be readily solved by the method of suc- 
cessive over-relaxation with -ff(^E') determined by iteratively solving equations (1) 
- (3) (see [12]). Figure 1 shows a sample solution exhibiting a Poynting jet. The 
solution consists of a region near the axis which is magnetically collimated by the 
toroidal field and a region far from the axis, on the outer radial boundary, which 
is anti- collimated in the sense that it is pushed against the outer boundary. The 
field lines returning to the disk at r > tq are anti-coUimated by the pressure of the 
toroidal magnetic field. 

Most of the twist A0 of a field line of a Poynting jet occurs along the jet from 
z = to Zm- Because —r^d(f)/H{'^) = dz/B^, we have 

A0(^) (^ot)-^(^/^o) Z, 



■'m 



(4) 



where r^Bz{r,z) is evaluated on the straight part of the jet at r = r(\E'). In 
the core of the jet \E' <C "^q, we have ~ 3^^^{'^q/'^)^^'^, and in this region we 
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FIGURE 1. Poloidal field lines for Poynting jet case for twist ujot = 1.84 rad. and 
{—H)max = l-13^'o/?'o with ^' =const contours measured in units of 5*0 which is the maxi- 
mum value of The outer boundaries at r = Rm and z = Zm are perfectly conducting and 
correspond to an external plasma. This external plasma will expand in response to the magnetic 
pressure of the jet so that Rm and Zm will increase with time. The dashed contour is the separa- 
trix with the X-point indicated. The initial dipole-like poloidal magnetic field is characterized by 
the flux function on the disk surface (r, z = 0) = {rl /2^/'^)r'^ K[l - {r / r^-lf /SI]/ {rl/2 + r'^f/'^ , 
where tq is the radius of the O— point of the field in the disk. In this and subsequent plots 5* is 
measured in units of = rlK/^^/"^. Note that 5^(0) « 10.4*o/?-o- 




can take * = C'^Qir/roY, and H = -/C(^oAo)(^/^o)*, where C, g, /C, and s 
are dimensionless constants. Equation (1) for the straight part of the jet imphes 
q = 1/(1 - s) and ^^(i-*-) = s(l - s)2/CV(l - 2s). Thus we find s = 1/4 so that 
q = 4/3, C = [9/32]2/3/C4/3, and /C = 3i/M(roCJot/^m). 

In order to have a Poynting jet, we find that /C must be larger than ?a 0.5. For 
the case of uniform expansion of the top boundary, = Vzt, this condition is the 
same as < 9.2(roc<Jo). For Figure (1), /C ^ 0.844. The field components in the 
straight part of the jet are 



These dependences agree approximately with those found in numerical simulations 

of Poynting jets [9]. On the disk, \E' ^ 3'^^'^^o{r/ro)'^ for r < ro/3^/'^. Using this 
and the formula for \E'(r) gives the relation between the radius of a field line in 
the disk, denoted r^, and its radius in the jet, r/r-Q = 6.5{rd/ro)^^'^]C~^. Thus the 
power law for ^ is applicable for vi < r < r2, where ri = 6.5ro(rj/ro)^/^//C and 
r2 = 1.9ro//C, with the inner radius of the disk. The outer edge of the Poynting 
jet has a transition layer where the axial field changes from Bz{r2) to zero while 
(minus) the toroidal field increases from —B^{r2) to {—H)max/i^2- Using equations 
(5), which are only approximate at r2, gives {—H)max ~ l-2/C\E'o/ro, which agrees 
approximately with our Grad-Shafranov solutions. 



Full, axisymmetric MHD simulations of the evolution of the coronal plasma of a 
Keplerian disk initially threaded by a dipole-like magnetic field are shown in Figure 
2. For these simulations the outer boundaries at r = Rm and z — Z^n were treated 
as free boundaries following the methods of [14]. These simulations established 
that a quasi- stationary collimated Poynting jet arises from the inner part of the 
disk while a steady uncoUimated hydromagnetic outfiow arises in the outer part of 
the disk. 

Quasi-stationary Poynting jets from the two sides of the disk within tq give an 
energy outfiow per unit radius of the disk dEs/dr = rvK{—B^Bz)h, where the 
h subscript indicates evaluation at the top surface of the disk. This outfiow is 



~ rodEB/dr ~ VK{ro){'^o/rof ~ W^'\eTg/s)rl{l^[B,{0)/6kGf where row is in 



units of lO^^cm, and Mg in units of IO^Mq. This formula agrees approximately 
with the values derived from the simulations. The jets also give an outfiow of 
angular momentum from the disk which causes disk accretion (without viscosity) 
at the rate Msir) = -2T:Evr = -2{r^ /vK){B^B^)h ~ 2^'g/[r^vx(ro)], where E is 
the disk's surface mass density [14]. The Poynting jet has a net axial momentum 
flux = (1/4) / rdr{Bl - Bfj ~ 0.5(*o/'^o)^ which acts to drive the jet outward 
through an external medium. Further, the Poynting jet generates toroidal magnetic 
flux at the rate $t ~ -12[vx(ro)/ro]^o- 





MHD SIMULATIONS OF POYNTING JETS 



For long time-scales, the Poynting jet is of course time-dependent due to the 
angular momentum it extracts from the inner disk (r < tq). This loss of angular 
momentum leads to a "global magnetic instability" and collapse of the inner disk 
[13]. An approximate model of this collapse can be made if the inner disk mass Md 
is concentrated near the O— point radius ro(t), if the field line slippage through the 
disk is negligible [13], \E'o =const, and if {—rB^)max ~ ^o/foif) (as found here). 
Then, Madro/dt = -2^g(G'Mro)"^/^. If U denotes the time at which ro(ti) = 
(the inner radius of the disk), then ro{t) = rj[l — (t — U) /tcou]'^^^ , for t < tj, where 

tcoii = V GM MdT^^'^/ (3\E'o) is the time-scale for the collapse of the inner disk. 
(Note that the time-scale for tq to decrease by a factor of 2 is ~ tj(r'o/rj)^/^ tj 
for vq ^ Tj.) The power output to the Poynting jets is E(t) = {2/3){AEtot/tcoii)[^ — 
(t - ti)/tcou]~^^^, where AEtot = GMMd/^n is the total energy of the outburst. 
Roughly, tcou ~ 2 dayM|(Md/M0)(6 x 1032QcJn2^^^^)2 ^ Schwarzschild black 
hole, where validity of the analysis requires tcoii 3> U. Such outbursts may explain 
the flares of active galactic nuclei blazar sources [16,17] and the one-time outbursts 
of gamma ray burst sources [18]. 




FIGURE 2. Time evolution of dipole-like field threading the disk from the initial configuration 
t = (bottom panels) to the final quasi-stationary state t = 1.2tout, where tout is the rotation 
period of the disk at the outer radius Rm of the simulation region from [9] . The left-hand panels 
show the poloidal field lines which arc the same as '^{r,z) =const lines; ^' is normalized by 
(the maximum value of ^'(r, z)) and the spacing between lines is 0.1. The middle panels show the 
poloidal velocity vectors Vp. The right-hand panels show the constant lines of —rB^{r, z) > in 
units of '^o/ro and the spacing between lines is 0.1. 



CONCLUSIONS 



Recent MHD simulation studies support the idea that an ordered magnetic field 
of an accretion disk can give powerful outflows or jets of matter, energy, and angular 
momentum. Most of the studies have been in the hydromagnetic regime and find 
asymptotic flow speeds of the order of the maximum Keplerian velocity of the disk, 
VK'i- These flows are clearly relevant to the jets from protostellar systems which 
have flow speeds of the order of Vki- In contrast, observed VLBI jets in quasars and 
active galaxies point to bulk Lorentz factors F ~ 10 - much larger than the disk 
Lorentz factor. In the jets of gamma ray burst sources, F ~ 100. The large Lorentz 
factors as well as the small Faraday rotation measures point to the fact that these 
jets are in the Poynting flux regime. These jets may involve energy extraction from 
a rotating black hole [19, 20]. 
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